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1.  INTRODUCTION 


This  work  has  been  a  theoretical  study  of  the  process  of  acoustic  scattering 
from  small  objects  at  the  water-sediment  interface  or  buried  in  the  sediment.  An 
object  is  regarded  as  "small"  if  its  dimensions  are  small  compared  to  the 
acoustic  wavelengths  in  the  frequency  band  of  interest.  At  normal  incidence, 
acoustic  scattering  from  these  objects  is  a  fairly  simple  problem,  which  has 
been  analyzed  in  connection  with  depth  sounder  applications.  Our  study  is 
aimed  at  the  shallow  grazing  angie  case,  which  is  particularly  relevant  to 
detection  applications.  In  particular,  we  would  like  to  obtain  estimates  of  the 
degradation  in  the  acoustic  backscatter  of  the  buried  objects.  We  will  consider 
acoustic  frequencies  less  than  1  kHz  such  that  the  acoustic  wavelength  will  be 
large  compared  to  objects  of  about  a  meter  in  size. 

This  report  is  divided  into  five  main  sections:  In  Section  2,  the  theoretical 
tools  will  be  described.  Both  the  physical  models  and  the  systems  parameters 
will  be  described.  In  Section  3,  the  application  of  reciprocity  to  estimate  the 
backscattered  signal  is  described.  In  Section  4,  the  resulting  theoretical  model 
is  applied  to  specific  cases.  The  results  and  conclusions  are  summarized  in 
Section  5. 
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2.  THEORY 


The  theoretical  considerations  may  be  divided  into  two  main 
complementary  parts,  physical  processes  and  systems  analysis.  The  former 
concerns  the  sound  propagation  and  scattering  models  and  the  latter  concerns 
the  system  parameters,  particularly  the  degradation  in  target  echo. 

2.1  PHYSICAL  PROCESSES 

We  require  a  theoretical  model  that  can  accurately  predict  the  propagation 
of  sound  waves  through  the  water,  through  the  water/sediment  interface,  the 
interaction  of  the  sound  waves  with  a  small  object,  and  the  backscattered  signal 
from  the  object  back  to  the  point  of  transmission.  Analytical  solutions  are  only 
practical  for  plane  wave  problems  with  one  interface.  Clearly,  this  will  not  be 
adequate  for  studying  the  backscatter  of  sound  from  small  objects,  particularly 
in  cases  where  the  object  is  near  or  at  the  interface,  or  where  the  medium  is 
multilayered.  Therefore  it  was  decided  to  approach  the  problem  through 
numerical  computer  simulations. 

Objects  may  become  buried  in  sediment  through  a  variety  of  mechanisms, 
including  sedimentation  of  waterborne  material  or  scouring  of  the  surrounding 
material;  or,  the  object  may  simply  sink  into  a  soft  sediment  under  its  own 
weight.  This  type  of  burial  implies  either  sand  or  silt  sediments,  and  thus  the 
study  will  be  limited  to  these  two  types  of  sediments.  The  detection  of  objects 
that  are  embedded  in  solid  rock  sediment  will  not  be  considered. 

Sand  and  silt  sediments  are  porous  solid  media.  The  theory  of  sound 
propagation  in  porous  media  was  first  developed  by  Biot  The  porous  solid 
differs  from  the  viscoelastic  solid  in  that  it  can  support  two  longitudinal  waves,  a 
fast  wave  and  a  slow  wave.  The  fast  wave  corresponds  to  the  acoustic  wave  in 
a  viscoelastic  solid,  but  the  slow  wave  has  no  counterpart  in  a  viscoelastic  solid. 
Stern®  has  developed  a  numerical  model  for  a  horizontally  stratified  porous 
medium,  but  it  only  deals  with  plane  waves;  and  Norris7  has  published  a 
theoretical  model  for  radiation  and  scattering  of  sound  in  a  homogeneous 
porous  medium.  Thus,  it  appears  that  most  of  the  compon/'-'ts  for  producing  a 
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propagation  model  in  porous  media  are  available,  but  they  have  not  yet  been 
assembled  in  a  comprehensive  package. 

The  principal  difference  between  a  porous  solid  and  a  viscoelastic  one  is 
the  presence  of  the  slow  wave.  Stern6  has  shown  that  the  slow  wave  velocity  is 
frequency  dependent.  From  his  results,  it  appears  that  for  frequencies  below 
1  kHz,  the  slow  wave  velocity  is  too  low  to  have  any  significant  effect  on  the 
sound  field,  and  thus  the  sediment  may  be  approximated  as  a  viscoelastic  solid. 
Therefore,  it  is  not  imperative  that  we  use  the  Biot  model  for  porous  media. 
Since  mathematical  models  for  viscoelastic  solids  already  exist  in  very 
comprehensive  packages,  we  have  chosen  to  approximate  the  sediment  as  a 
viscoelastic  solid. 

For  computer  simulation  of  sound  propagation  in  viscoelastic  media,  the 
finite  difference8  and  finite  element0  methods  are  the  most  general  methods  but 
their  computation  resource  requirements  are  prohibitive.  The  parabolic 
equation  method10  is  inadequate  because  it  ignores  shear  stress  and  strain. 
Ray  tracing  methods11  are  also  inadequate  because  they  are  only  high 
frequency  approximations  and  are  expected  to  be  particularly  inaccurate  in  the 
vicinity  of  the  interfaces. 

Of  the  computer  models  that  are  currently  available,  the  seismic  and 
acoustic  fast-field-program  for  range  independent  environments  (SAFARI) 
model,  developed  at  the  S^OLANT  Research  Centre  by  Schmidt  and 
Jensen,1^-18  is  the  most  complete.  It  gives  a  full  wave  solution  for  the  field 
generated  by  a  point  source,  in  a  horizontally  stratified  viscoelastic  medium. 
This  is  done  by  decomposing  the  acoustic  field  from  a  point  source  into  a  series 
of  plane  waves,  and  then  solving  the  field  solutions  for  all  plane  wave 
components  and  superimposing  them  to  give  the  total  field.  The  algorithm  is  an 
integral  part  of  the  signal  processing  package  that  can  simulate  the  propagation 
of  an  arbitrary  signal  projected  from  an  arbitrary  array  of  point  sources. 
Therefore  the  SAFARI  model  will  be  our  main  model  fo<  quantitative  analyses. 

For  a  qualitative  look  at  the  acoustic  pressure  field,  the  complex  ray 
model14  is  very  efficient.  It  is  particularly  efficient  when  applied  to  Gaussian 
beams.  Results  from  the  complex  ray  model  have  been  successfully  compared 
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with  the  results  of  the  full  wave  methods  in  bench  mark  tests.  In  comparisons 
with  the  SAFARI  model,  significant  quantitative  difterences  were  found  but, 
qualitatively,  their  results  were  very  similar.  That  is,  the  two  models  showed 
similar  acoustic  pressure  amplitude  distributions,  but  there  were  differences  in 
the  absolute  levels.  The  complex  ray  model  only  takes  a  small  fraction  of  the 
computation  time  that  SAFARI  takes.  The  computation  time  requirement  of  the 
complex  ray  model  is  relatively  insensitive  to  the  angle  of  incidence,  while  that 
of  the  SAFARI  model  increases  rapidly  with  grazing  angle.  In  fact,  it  is  simply 
not  practical  to  use  the  SAFARI  for  near-normal  incidence  angles.  Therefore, 
the  complex  ray  model  will  be  used  as  a  qualitative  support. 

The  main  deficiency  in  the  SAFARI  model  is  the  constraint  of  a  horizontally 
stratified  medium.  This  means  that  SAFARI  cannot  directly  model  the  signal 
scattered  from  a  small  object.  It  may  be  used  to  predict  the  propagation  loss  to 
and  from  the  object,  but  the  scattering  strength  of  the  object  itself  must  be 
computed  separately. 

There  are  very  few  experimental  verifications  of  the  SAFARI  model.  An 
experiment  has  been  conducted  recently  by  Satkowiak,15  Naval  Coastal 
Systems  Center  (NCSC),  in  cooperation  with  the  Admiralty  Research 
Establishment  (ARE),  U.K.,  at  the  ARE  facility  in  Portland,  U.K.  Experiments 
conducted  by  Applied  Research  Laboratories,  The  University  of  Texas  at  Austin 
(ARL:UT)  in  cooperation  with  other  Navy  laboratories,  off  Kings  Bay,  Georgia, 
and  in  cooperation  with  NCSC,  off  Panama  City,  gave  mixed  results.  These 
were  high  frequency  acoustic  experiments,  in  the  region  of  10-100  kHz,  where 
the  slow  wave  velocity  is  such  that  significant  coupling  may  occur.  For  the 
frequency  range  of  interest,  that  is,  less  than  1  kHz,  as  far  as  we  know,  the 
viscoelastic  solid  approximation  is  valid. 

2.2  SYSTEMS  ANALYSIS 

Let  us  start  with  the  sonar  equation.  As  illustrated  in  Fig.  2.1,  the  signal  S{0 
returned  by  a  object  of  target  strength  Ts  at  a  range  r  in  an  unbounded  body  of 
water  is  given  by 


(a)  DETECTION  IN  WATER 


(b)  DETECTION  IN  SEDIMENT 


FIGURE  2.1 

SYSTEM  PARAMETERS  FOR  TARGET  DETECTION 

ARL  -  LIT 

AS-89-9 
NPC - RFO 
1  -  20  -  89 
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Sto  =  S0  +  Ts  -  l_2ro 


(2-1) 


where  S0  is  the  source  level  and  L2ro  is  the  two-way  propagation  loss  in  water 
at  the  range  r.  For  both  source  and  target  in  a  homogeneous  medium,  the  two- 
way  propagation  loss  is  simply  twice  the  one  way  loss  Lro, 

L2ro  =  2  Lro  .  (2-2) 

Let  us  use  the  above  equations  as  our  standard  against  which  the  return  from 
small  targets  at  the  sediment  interface  or  buried  within  the  sediment  will  be 
compared.  With  reference  to  Fig.  2.1,  the  return  St  from  the  same  target  at  the 
interface  or  within  the  sediment  is  given  by 

St  =  S0  +  Ts  -  L2r  ,  (2-3) 

where  L2r  is  the  effective  two-way  propagation  loss. 

Our  objective  is  to  estimate  the  change  in  signal  level  due  to  the  presence 
of  the  sediment,  which  we  will  call  the  "sediment  factor"  Fs,  defined  as 

Fs  =  St  -  Sto 

=  L2r  '  1-2 ro  .  (2-4) 

It  is  a  quantitative  measure  of  the  signal  degradation  caused  by  the  presence  of 
the  sediment. 

Let  us  examine  L2r,  the  main  constituent  of  the  sediment  factor,  more 
closely.  It  is  implied  in  the  above  equations  that  the  target  strength  Ts  is  the 
target  strength  as  measured  in  water.  If  the  object  were  placed  in  the  sediment, 
then  its  target  strength  will  undergo  a  change  Lz.  It  is  also  possible  that  the 
pressure  propagation  loss  along  the  outward  bound  path  Lr+  may  not  be  equal 
to  that  of  the  return  path  Lr..  The  effective  two-way  propagation  loss  may  be 
defined  as  the  sum  of  these  three  components, 

*-2r  =  *-r+  +  Lr-  +  Lz  (2-5) 
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Due  to  the  complicated  nature  of  the  propagation  process,  we  will  rely  on 
computer  simulations  to  give  the  quantitative  estimates  of  Lr+.  Assuming  that 
return  paths  are  simply  the  reverse  of  the  outward  bound  paths,  it  is  possible  to 
estimate  the  propagation  loss  difference  Lr.  -  Lr+  by  reciprocity  and  the  target 
loss  Lz  from  the  sediment  properties. 
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3.  THE  BACKSCATTERED  SIGNAL 


3.1  TARGET  STRENGTH 

The  target  strength  of  small  objects  are  of  immediate  interest.  Target 
strength  Ts  in  decibels  is  defined  in  terms  of  the  acoustic  cross  section  o  by 

Ts  =  10log10(o/4jt)  .  (3-1) 

For  a  solid  sphere  that  is  small  compared  to  the  wavelength,  the  acoustic  cross 
section  derived  by  Rayleigh16  and  later  by  Anderson1 7  is  given  by 


^sphere  =  sP  (ka)^  R  ,  (3-2) 

where  k  is  the  acoustic  wave  number  in  the  surrounding  medium  and  a  is  the 
radius;  R  is  a  constant  dependent  on  material  properties, 

R  - 10  -gh2)/(3gh2)  +  (1  -g)/(1  +2g)]2  ,  (3-3) 

where  h  and  g  are  the  sound  velocity  and  density  ratios  between  the  sphere 
and  the  surrounding  medium. 

Let  us  examine  a  few  asymptotic  cases:  consider  the  cases  where  the 
values  of  g  (velocity  ratio)  and  h  (density  ratio)  are  either  0,  1 ,  or  infinite.  The 
resulting  values  of  R  for  all  reasonable  combinations  is  given  in  the  table  below. 
Certain  combinations,  such  as  h=°°  and  g=0,  are  not  reasonable  since  they  do 
not  yield  deterministic  results  and  do  not  approximate  any  object  constructed  of 
real  materials. 
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TABLE  3.1 

ASYMPTOTIC  VALUES  OF  R 


g 

D 

R 

Remarks 

0 

0 

oo 

Trapped  gas  bubbles 

i 

4 

0 

Undetectable  object 

i 

OO 

0.33 

Very  rigid  object  of  same  density  as  medium 

oo 

1 

0.69 

Very  dense  object  of  same  sound  velocity  as 
the  medium 

oo 

OO 

0.69 

Very  rigid  and  very  dense  object 

It  is  interesting  to  note  that  a  low  rigidity  and  low  density  object  will  have  a 
very  large  value  of  R,  indicating  that  trapped  gas  bubbles  may  be  very  strong 
sound  scatterers.  An  object  with  the  same  density  and  sound  speed  as  the 
medium  will  not  be  detectable.  In  water,  silt,  or  sand,  there  is  a  large  class  of 
man-made  objects  that  can  be  considered  as  being  either  very  much  more 
dense  or  very  much  more  rigid  than  the  medium.  It  can  be  seen  that  such 
targets  will  have  R  values  between  0.33  and  0.69. 

Rewriting  Eq.  (3-2)  in  terms  of  volume,  we  get 

asphere  =  (9/4rc)  R  ,  (3-4) 

where  V  is  the  volume  of  the  sphere. 

It  is  expected  a  similar  expression  would  be  applicable  for  objects  that  are 
very  small  compared  to  the  wavelength.  It  is  known,  for  example,  that  for  any 
dense  and  rigid  elongated  solid  of  revolution,1®  the  acoustic  cross  section  in 
the  axial  direction  is  independent  of  shape,  given  by 

arigid  dense  elongated  =  filn)  ^  -  (3'®) 

where  V  is  the  volume  of  the  object. 
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It  is  seen  that  Eqs.  (3-4)  and  (3-5)  are  of  the  same  form,  differing  only  in  the 
values  of  the  constant  multiplier. 

Therefore,  the  acoustic  backscattering  cross  section  of  small  objects  can 
be  expected  to  vary  as 

o  =  C  V2  k4  ,  (3-6) 

where  the  constant  C  is  weakly  dependent  on  the  shape  and  orientation  of  the 
object. 

It  is  seen  that  target  strength  is  proportional  to  the  fourth  power  of  the 
speed  of  sound  in  the  surrounding  medium.  Therefore,  the  change  in  target 
strength  incurred  by  moving  a  target  from  one  medium  into  another  is  simply 
given  by  the  fourth  power  of  the  sound  velocity  ratio, 

Lz  =  -  40  log10(rc)  ,  (3-7) 

where  rc  is  the  sound  velocity  in  the  second  medium  divided  by  the  sound 
velocity  in  the  first  medium. 

3.2  RECIPROCITY 

All  of  the  above  models  assume  that  the  medium  is  linear.  Therefore  the 
signal  path  between  source  and  receiver  may  be  considered  as  a  linear 
network.  From  linear  network  theory,  it  is  known  that  the  forward  and  reverse 
paths  must  satisfy  the  reciprocity  theorem.  In  this  context,  the  reciprocity 
theorem  may  be  interpreted  as  follows:  If  a  source  at  a  position  A,  defined  in 
terms  of  volume  displacement,  produces  a  certain  level  of  pressure  at  a 
position  B,  then  a  source  of  equal  strength  at  B  will  produce  the  same  pressure 
at  A. 


Both  the  SAFARI  and  complex  ray  models  estimate  the  pressure 
propagation  at  any  point  B  produced  by  a  unit  pressure  source  at  A,  where  the 
source  level  is  defined  in  terms  of  pressure  at  a  reference  radial  distance  of  1  m. 
Therefore  the  pressures  at  A  and  B  produced  by  unit  sources  at  B  and  A, 
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respectively,  are  expected  to  differ  by  a  factor  equal  to  the  ratio  of  the  acoustic 
impedances  at  the  two  points.  Therefore  the  propagation  loss  difference  is 
given  by 


Lr_-Lr+  =  20log10(rz)  ,  (3-8) 

where  rz  is  the  local  acoustic  impedance  at  the  target  position  divided  by  that  at 
the  source  position. 

The  acoustic  impedance  in  an  homogeneous  medium  is  equal  to  the 
product  of  the  density  and  the  sound  speed.  This  is  also  approximately  true  in  a 
layered  medium  for  points  several  wavelengths  away  from  the  discontinuities. 
For  points  near  or  on  a  boundary  between  two  layers,  the  impedance  is  not 
easily  determined,  but  it  is  expected  to  have  an  intermediate  value.  Thus,  let  us 
approximate  the  effective  impedance  Z,  at  a  point  in  a  medium  with  impedance 
Zi  at  a  distance  d  from  the  interface  with  a  medium  of  impedance  Z2,  by 

Z  =  (Z^^l  +Z2)/(e'd/xl  +  1)  .  (3-9) 

where  X.-|  is  the  wavelength  in  the  medium  with  impedance  Zi . 
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4.  COMPUTER  SIMULATIONS 
4.1  SEDIMENT  PROPERTIES 

As  stated  earlier,  the  study  is  limited  to  sediments  that  may  be  described  as 
either  silt  or  sand,  or  some  combination  thereof.  For  practical  purposes,  let  us 
restrict  the  study  to  the  first  100  m  of  the  sediment.  Over  such  a  short  depth 
interval,  there  is  usually  very  little  change  in  the  sediment  properties,  and 
therefore  they  may  be  regarded  as  constant.  A  set  of  typical  values1 9'23  of  the 
properties  of  silt  and  sand  in  the  upper  100  m  is  shown  in  Table  4.1.  These 
values  will  be  used  in  the  numerical  models  to  compute  the  sediment  factor  of 
signals  backscattered  by  a  small  object.  In  these  simulations,  the  sonar  is 
positioned  at  a  height  of  100  m  above  the  sediment  interface.  The  sound 
velocity  and  density  of  the  water  are  assumed  to  be  1540  m/s  and  1000  kg/m3, 
respectively. 


TABLE  4.1 

TYPICAL  SEDIMENT  PROPERTIES 


Silt  sediment 


Density 

1400 

kg/m 

Sound  velocity 

1500 

m/s 

Sound  velocity  gradient 

1 

s*1 

Sound  absorption  coefficient 

0.5 

dB/X 

Shear  wave  velocity 

0 

m/s 

Shear  absorption  coefficient 

- 

dB/X 

Sand  sediment 

Density 

2000 

kg/m 

Sound  velocity 

1725 

m/s 

Sound  velocity  gradient 

0 

s'1 

Sound  absorption  coefficient 

0.5 

dB/X 

Shear  wave  velocity 

55 

m/s 

Shear  absorption  coefficient 

3.44 

dB/X 
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4.2  ACOUSTIC  PRESSURE  FIELD 


In  a  water  logged  silt  sediment,  the  sound  velocity  is  almost  identical  to  that 
cf  the  water.  In  gassy  silt  sediments,  the  sound  velocities  are  often  less  than 
that  of  water,  in  these  sediments,  sound  propagates  directly  into  the  sediment 
by  refraction  and  with  very  little  loss.  There  is  a  critical  angle  in  the  sediment 
which  forces  the  refracted  waves  to  emerge  at  incidence  angles  less  than  or 
equal  to  the  critical  angle.  Beyond  a  corresponding  critical  range,  this  effect 
becomes  dominant.  There  is  the  possibility  of  interference  effects  between  the 
sound  waves  from  the  two  sides  of  the  critical  range,  as  illustrated  in  Fig.  4.1 . 
Thus,  the  acoustic  field  may  be  considered  in  terms  of  three  zones,  the 
refraction,  critical  refraction,  and  interference  zones. 

In  a  compacted  sand  sediment,  the  sound  velocity  is  significantly  higher 
than  that  of  water.  Acoustic  energy  may  be  introduced  into  the  sediment  by  two 
mechanisms,  refraction  and  tunneling.  Near  normal  incidence,  sound  waves 
are  refracted  into  the  sediment  in  the  usual  way.  At  the  critical  angle,  the  sound 
energy  is  refracted  in  a  direction  that  is  almost  parallel  to  the  interface.  For 
incidence  angles  greater  than  the  critical  angle,  there  is  total  reflection.  Thus, 
beyond  a  corresponding  critical  range,  sound  may  not  enter  the  sediment  by 
refraction.  Instead,  it  enters  as  an  evanescent  wave  that  tunnels  into  the 
sediment  to  a  depth  of  a  few  wavelengths.  Since  the  refracted  and  evanescent 
waves  travel  at  different  velocities  they  are  expected  to  produce  an  interference 
pattern.  Thus,  the  acoustic  field  may  be  considered  in  terms  of  three  zones,  the 
refraction,  the  evanescent,  and  the  interference  zones,  as  shown  in  Fig.  4.2. 

A  numerical  example  of  the  sound  pressure  contours  for  a  sandy  sediment, 
computed  by  the  complex  ray  model,  is  shown  in  Fig.  4.3.  It  shows  the 
in-sediment  constant  pressure  amplitude  contours  produced  by  a  unit  point 
source  100  m  above  the  sediment  operating  at  100  Hz.  The  critical  angle  is 
calculated  to  be  63.2°;  the  critical  ray  intersects  the  interface  at  a  horizontal 
distance  of  198  m.  It  is  seen  that  for  ranges  greater  than  about  200  m,  the 
sound  field  appears  to  be  divisible  into  three  distinct  regions.  Nearest  the 
interface,  there  is  a  region  where  the  contours  of  constant  sound  pressure  are 
nearly  horizontal  and  very  tightly  packed,  indicative  of  the  rapid  decay  rate  of 
the  evanescent  pressure  with  depth.  Farthest  from  the  interface,  there  is  a 
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FIGURE  4.1 

SOUND  PROPAGATION  FROM  WATER  INTO  A  SILTY  SEDIMENT 


ARL-UT 
AS-89-10 
NPC - RFO 
1-20-89 
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FIGURE  4.2 

SOUND  PROPAGATION  FROM  WATER  INTO  A  SANDY  SEDIMENT 


ARL - UT 
AS-89-1 1 
NPC  -  RFO 
1-20-89 


region  of  near-vertical  contours  that  are  widely  spaced.  In  between,  there  is  a 
region  in  which  the  contours  change  direction  and  form  closed  loops  around 
periodic  nulls.  These  three  regions  correspond  to  the  refraction,  evanescent, 
and  interference  zones  indicated  in  Fig.  4.2.  The  periodic  nulls  are  the  result  of 
interference  between  the  evanescent  and  refracted  waves.  The  spacing 
between  the  nulls  is  equal  to  the  spatial  period  of  the  interference  pattern. 

Another  example  was  computed  for  an  acoustic  frequency  of  50  Hz,  as 
shown  in  Fig.  4.4.  It  shows  a  similar  pressure  distribution,  but  the  spatial  period 
of  the  interference  nulls  are  twice  that  of  the  100  Hz  case. 

The  above  pressure  contours  were  computed  using  the  complex  ray 
model,  which  assumes  a  fluid-fluid  medium.  The  complex  ray  model  is  a  very 
efficient  algorithm,  and  therefore  it  is  useful  for  obtaining  qualitative  estimates  of 
the  sound  field.  For  quantitative  analysis,  we  use  the  SAFARI  model,  which  is  a 
full  wave  solution  for  viscoelastic  solid  and  liquid  media.  From  the  contours  in 
Figs.  4.3  and  4.4,  it  is  seen  that  the  most  interesting  region  appears  to  be 
beyond  a  range  of  300  m  and  above  a  depth  of  50  m.  Therefore,  to  conserve 
computer  resources,  subsequent  computations  were  confined  to  this  region. 

The  pressure  fields,  as  calculated  by  the  SAFARI  model  for  a  point  source 
100  m  above  the  interface  projecting  a  100  Hz  acoustic  wave,  are  shown  in 
Figs.  4.5  and  4.6  for  the  silt  and  sand  sediments.  These  contours  should  be 
compared  with  the  pressure  contours  for  the  case  of  no  sediment  at  all,  shown 
in  Fig.  4.7. 

It  is  seen  that  the  propagation  loss  in  the  silty  sediment  is  not  significantly 
different  from  that  of  the  no-sediment  case  up  to  a  range  of  about  900  m. 
Beyond  this  range,  the  effects  of  the  critical  angle  in  the  sediment  become 
apparent.  There  is  an  interference  zone  corresponding  to  the  one  illustrated  in 
Fig.  4.1,  in  which  there  are  significant  pressure  minima  and  maxima. 

The  pressure  in  the  sand  sediment  computed  by  the  SAFARI  model 
appears  to  be  similar  in  character  to  that  produced  by  the  complex  ray  model 
shown  in  Fig.  4.3,  and  the  illustration  in  Fig.  4.2,  showing  the  evanescent, 
refraction,  and  interference  zones.  At  a  more  detailed  level,  the  SAFARI  model 
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appears  to  give  a  more  complicated  pressure  field  than  the  complex  ray  model. 
This  is  not  altogether  surprising  since  there  are  significant  differences  between 
the  two  models;  the  main  shortcoming  of  the  complex  ray  model  is  that  it  is  only 
a  fluid-fluid  model.  A  discussion  of  the  differences  and  their  reasons  are 
beyond  the  scope  of  this  report. 

4.3  SEDIMENT  FACTOR 

Let  us  now  consider  a  sonar  with  a  widebeam  projector  that  could  be 
approximated  by  a  point  source,  and  a  similar  widebeam  hydrophone. 

Since  a  unit  source  was  used,  the  inverse  of  the  pressure  contours  in 
Figs.  4.5  and  4.6  may  be  directly  equated  to  the  forward  propagation  loss  Lr+. 
Using  Eqs.  (3-7)  and  (3-8),  the  loss  terms  Lr_  and  Lz  were  estimated. 
Substituting  into  Eq.  (2-5),  the  total  propagation  loss  L2r  was  estimated.  Finally, 
substituting  into  Eq.  (2-4),  the  sediment  factor  Fs  was  computed.  The  results  for 
the  silt  and  sand  sediments  are  shown  in  Figs.  4.8  and  4.9. 

In  Fig.  4.8,  it  is  seen  that  the  sediment  factor  in  the  silty  sediment  is  quite 
small  for  ranges  up  to  1000  m.  Beyond  this  range,  the  interference  effects  due 
to  the  critical  angle  in  the  sediment  begin  to  significantly  distort  the  field, 
producing  a  patch  of  gain  greater  than  5  dB  at  a  range  of  1 100  m  and  a  depth  of 
25  m,  and  then  giving  way  to  increasing  losses  at  longer  ranges. 

The  sediment  factor  contours  in  the  sand  sediment,  shown  in  Fig.  4.9, 
follow  the  same  trends  as  the  one-way  propagation  loss  contours  shown  in 
Fig.  4.6.  The  refracted,  evanescent,  and  interference  zones  are  clearly  visible. 
At  100  Hz,  the  acoustic  wavelength  is  approximately  17  m.  In  the  evanescent 
zone,  it  is  seen  that  the  sediment  factor  drops  to  a  value  less  than  -30  dB  at  a 
depth  of  one  wavelength.  The  -30  dB  contour  appears  to  coincide  with  the 
boundary  between  the  distinct  zones.  Beyond  a  range  of  400  m,  the  -30  dB 
contour  marks  the  boundary  between  the  evanescent  wave  zone  and  the 
interference  zone.  Below  a  depth  of  approximately  one  wavelength,  it  marks  the 
division  between  the  refraction  zone  and  the  interference  zone. 
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The  sediment  factor  contours  may  be  plotted  in  a  dimensionless  space  that 
will  allow  them  to  be  applied  to  a  broader  range  of  operating  conditions  than  the 
specific  ones  modeled  here.  We  could  have  chosen  three  dimensionless 
parameters;  that  is,  target  depth,  horizontal  range,  and  sonar  height  in  acoustic 
wavelengths.  The  first  must  be  included  because  it  determines  the  rate  of 
decay  of  evanescent  waves.  The  remaining  two  are  not  completely 
independent;  it  is  expected  that  they  may  be  replaced  by  a  single  parameter 
such  as  the  incidence  or  grazing  angle  because  of  the  strong  influence  the 
critical  angle  has  on  the  acoustic  field.  We  chose  the  grazing  angle.  The 
grazing  angle,  however,  is  difficult  to  define  in  a  situation  where  ray  paths  are 
not  clearly  discernible.  In  the  absence  of  a  better  definition,  we  have  chosen  to 
define  the  grazing  angle  as  the  inverse  tangent  of  the  sonar  height  divided  by 
the  horizontal  range.  The  results  are  shown  in  Figs.  4.10  and  4.1 1.  Although, 
strictly  speaking,  these  plots  will  change  with  the  operating  frequency  and  sonar 
height  above  bottom,  the  dependence  is  expected  to  be  weak,  particularly 
outside  the  interference  zones.  Therefore,  these  plots  are  expected  to  be  a 
useful  repiesentation  of  the  behavior  of  small  buried  targets  in  general. 
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5.  RESULTS  AND  CONCLUSIONS 


We  have  looked  at  the  behavior  of  a  small  target  in  two  representative 
sediment  types,  a  soft  silty  sediment  whose  sound  velocity  is  slightly  less  than 
that  of  water  and  a  sandy  sediment  with  a  much  higher  sound  velocity.  The 
reduction  in  the  target  echo  due  to  the  presence  of  the  sediment  was  quantified 
as  a  "sediment  factor".  Specifically,  the  sediment  factor  was  examined  both 
qualitatively  and  quantitatively.  Qualitatively,  it  was  found  that  in  both  cases,  the 
sound  field  in  the  sediment  may  be  divided  into  distinct  zones;  in  the  silty 
sediment  case,  they  are  the  refraction,  critical,  and  interference  zones,  and  in 
the  sand  sediment  case,  they  are  the  refraction,  evanescent,  and  interference 
zones.  For  quantitative  results,  computer  simulations,  using  the  SAFARI  full 
wave  mathematical  model,  were  used  to  obtain  numerical  estimates  of  the 
sediment  factor  for  a  widebeam  sonar  at  a  height  above  bottom  of  1 00  m  and 
operating  at  100  Hz.  In  a  silty  sediment,  the  sediment  factor  was  near  0  dB  in 
the  refraction  zone;  beyond  the  refraction  zone,  the  sediment  factor  varied  about 
a  declining  mean  due  to  the  critical  angle  and  interference  effects.  In  the  sand 
sediment,  the  sediment  factor  in  the  refraction  zone  was  mainly  caused  by  the 
absorption  coefficient  of  the  sediment;  in  the  evanescent  wave  zone  the 
sediment  factor  declined  with  increasing  depth,  reaching  -30  dB  at  a  depth  of 
only  one  wavelength;  in  the  interference  zone,  the  sediment  factor  varied 
periodically,  due  to  the  interference  between  the  refracted  and  evanescent 
waves,  about  a  mean  value  of  around  -50  dB.  Dimensionless  plots  of  the 
sediment  factor  as  a  function  depth  in  wavelengths  and  grazing  angle  in 
degrees  were  also  obtained.  These  plots  are  believed  to  give  a  useful 
representation  of  the  behavior  of  small  buried  targets  in  general.  The  sediment 
factor  due  to  the  presence  of  the  sediment  was  found  to  be  more  complicated 
than  expected.  The  results  obtained  in  this  report  are  based  entirely  on 
computer  models  and  they  remain  to  be  verified  by  an  experiment. 
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San  Diego,  CA  92152 

20 
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